Alternative determinism principle for topological analysis of chaos 
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The topological analysis of chaos based on a knot-theoretic characterization of unstable periodic 
orbits has proved a powerful method, however knot theory can only be applied to three-dimensional 
systems. Still, the core principles upon which this approach is built, determinism and continuity, 
apply in any dimension. We propose an alternative framework in which these principles are en- 
forced on triangulated surfaces rather than curves and show that in dimension three our approach 
numerically predicts the correct topological entropies for periodic orbits of the horseshoe map. 
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Chaotic behavior results from the interplay of two ge- 
ometrical processes in state space: stretching separates 
neighboring trajectories while squeezing maintains the 
flow within a bounded region A topological anal- 

ysis has been developed to classify the ways in which 
stretchingand squeezing can organize a chaotic attrac- 
tor Q, H 3 • It relies on a theorem stating that unsta- 
ble periodic orbits (UPO) of a chaotic three-dimensional 
(3D) flow can be projected onto a 2D branched manifold 
(a template) without modifying their knot invariants 
In this method, UPO extracted from an experimental 
time series are characterized by the simplest template 
compatible with their topological invariants 

However this approach can only be applied to 3D at- 
tractors: in higher dimensions, all knots can be deformed 
into each other. Although other topological methods are 
applicable to higher dimensions [(jljj, extending template 
analysis is still desirable because it provides a different 
information. A first step to overcome the 3D limitation is 
to recognize that knot theory is not a necessary ingredi- 
ent but simply a convenient way to study how two funda- 
mental properties, determinism and continuity, constrain 
trajectories in phase space. It is because two trajectories 
cannot intersect that the knot type of a 3D periodic orbit 
is well defined and is not modified as the orbit is deformed 
under control parameter variation. 

In this paper, we note that a dimension-independant 
formulation of determinism is orientation preservation 
and propose an approach where it is enforced on a repre- 
sentation of the dynamics in a triangulation of periodic 
points. In dimension three, an explicit formalism is easily 
constructed, and we find that it numerically predicts the 
correct entropies for periodic orbits of the horsehoe map. 
The entropy of a periodic orbit is an invariant defined as 
the minimal topological entropy @ of a flow containing 
this orbit E| ; a positive-entropy orbit is a powerful 
indicator of chaos 0, ITTl ll2T ] . This result suggests that a 
key ingredient for constructing a knotless template anal- 
ysis has been captured, although a proof of validity and 
an explicit higher-dimensional extension are still lacking. 

We now detail our approach. A first step is to replace 
the requirement of non-intersecting curves by a geomet- 



rical problem that adapts naturally to phase spaces of 
any dimension. It has been suggested to exploit the 
rigid structure of invariant manifolds of UPO 0, IT^ . 
Here, we note that when a volume element V of a d- 
dimensional phase space is advected by a deterministic 
flow <& t , the image ^t(dV) of its boundary cannot dis- 
play self-intersections: at any time t, its interior and its 
exterior remain distinct, as with a droplet in a fluid flow. 
A technical formulation of this property is that volume 
orientation is preserved by the dynamics. For simplicity, 
we consider attractors embedded in W 1 x S 1 (e.g., forced 
systems), that can be sliced into n-dimensional Poincare 
sections parameterized by <p G S 1 . Determinism then im- 
poses that boundaries of n-dimensional volume elements 
of Poincare sections retain their orientation (Fig. ^| . 




FIG. 1: Under the action of the flow, volume elements of 
Poincare sections and their boundaries are stretched and 
squeezed but retain their orientation, as illustrated here for 
2D section planes. 

Template analysis must be applicable to UPO ex- 
tracted from experimental signals, and thus can only rely 
on the phase space trajectory of a period-p orbit. Thus, 
we represent the dynamics in a triangulated space whose 
nodes are p periodic points Pi in a Poincare section, with 
Pi+i = F(Pi), F being the return map. In this space, 
points Pi are 0-cells, line segments (Pi, Pj) = (ij) joining 
two points are 1-cells, triangles (Pi,Pj,Pk) = (ijk) are 
2-cells, etc (Fig. EK) ■ Similar concepts have been used 
in |14| to analyze the static structure of an attractor, 
but we focus here on the dynamics. We denote by Sm 
the set of collections of contiguous m-cells, which are the 



analogs of m-dimensional surfaces in the original phase 
space. As Poincare sections are swept, periodic points 
move in the section plane and so do the m-cells attached 
to them (Fig. [2b) . The dynamics induced in S m should 
reflect that of m-dimensional phase-space surfaces under 
action of the chaotic flow, and in particular be organized 
by the same stretching and squeezing mechanisms. 




FIG. 2: (a) Triangulated space based on periodic points Pi 
in a 3D Poincare section. The 2-cell (ikl) is shaded, (b) The 
flow induces a mapping of this triangulated space into itself, 



as suggested here for a period-5 orbit embedded in 



s 1 . 



A dynamics in the triangulated space is specified by 
maps F m : S rn — > S m acting on collections of contiguous 
m-cells. Since the original return map F sends nodes to 
nodes but not facets to facets, the F m are not restrictions 
of F for m > 0. However we require them to mimic F 
in the following way: they should be invertible, satisfy 
determinism and result from a continuous deformation of 
facets, just as F is a continuous deformation of identity. 
The F m should also satisfy dF m (E) — P m _i(dE) where d 
is the boundary operator. As we see below, facets are not 
necessarily trivially advected between sections because 
degeneracies occur, at which action must be taken to 
preserve orientation. 

We now specialize to the 3D case. The volume element 
of a triangulated set of periodic points in a 2D Poincare 
section is a triangle (2-cell) based on three periodic points 
Pi, Pj, Pk- Let Pi(ip) be the position of Pi in section 



j i 

ip, with Pj(0) 
evolution of T 



Pi and P(2tt) = P i+l . The natural 
(Pi, P, , Pk) as ip increases is 



T(<p) = (P i (<p),P J (<p) i I\(<p)), 



(1) 



which would lead to a trivial induced return map F^ (T) = 
T(2n) = (Pj+i, Pj+i, Pfe+i) if Q was uniformly valid as 
a 2-cell. However, it is common that at some tp — tpo, 
one of the three points [say Pk(<p)] passes between the two 
others, thereby changing the orientation of the candidate 
2-cell T(<p) given by (Fig.[3J). As emphasized above, 
this is strictly forbidden by determinism, and we must 
thus modify the representation of the dynamics. It turns 
out that this problem has a simple solution. 

The degenerate triangle T(ipo) in Fig. [3] is like a flat- 
tened balloon whose boundary splits into two superim- 
posed sides with opposing outer normals. Determinism 



T(%) 




T(2n) 



FIG. 3: As three points move in the section plane when ip is 
increased, the triangle they form can change its orientation. 



is violated when these two sides cross each other so that 
interior and exterior, defined with respect to outer nor- 
mal, seem to be exchanged. However, the experimental 
data only constrain node motion, from which the facet 
dynamics is interpolated. To preserve determinism, we 
force the two opposing sides not to cross by swapping 
them at degeneracy, thereby canceling the inversion. 

This prescription is illustrated in Fig. 01 where the two 
opposing sides at triangle degeneracy are represented as 
a solid and a dashed line. The key point is that we con- 
struct the edge dynamics so that the left (solid line) and 
right (dashed line) sides remain at the left and right, 
respectively. Since the left (resp., right) side consists of 
itinerary (ik) + (kj) (resp., (ij)) before degeneracy and of 
itinerary (ij) (resp., (ik) + (kj)) after degeneracy, their 
relative position is preserved by applying the following 
dynamical rule in Si = S n —i at triangle inversion: 



(ik) + (kj) 



(ik) + (kj) 
(ij) 



(2a) 
(2b) 



These rules also apply to reverse paths (e.g., (ji) — » 
(jk) + (ki)). Note that dT = d(ijk) = (ij) + ((jk) + (ki)) 
is mapped by © to ({ik) + (kj)) + (ji) = d(ikj). The 
permutation compensates for triangle inversion so that 
orientation of dT, and hence determinism, is preserved. 




FIG. 4: A triangle is inverted as Pk passes between Pi and 
Pj. Identifying the solid (resp., dashed) paths in the initial 
and end configurations leads to substitution 

Itineraries visiting edges = (ij) in a given order 
are represented by words in a language A* over alphabet 
A = {eim}, and (J2J by an operator afj that in each word 
w replaces the letter ey by the string e^efej and eik&kj 
by e.y [hence (c*j ) 2 = 1]. For example, 



& ij &kl&li&ij&jl&li&ik£kj&ji • ■ ■ — ^■kl^li^ik^kj^jl^-li^-ij^jk^ki • • 
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The afj generate a non-trivial dynamics, as the image 
of an itinerary depends on how periodic points rotate 
around each other. This simple dynamics faithfully re- 
flects that of the flow around the periodic orbit, as we 
show by computing the entropy of the orbit. 

From the motion of periodic points Pi (ip) in the section 
plane as ip is swept, a list of I triangle inversions crf™ Jm 
is obtained, from which we build an induced return map 
that transforms a word w £ A* into another word w' as: 



m Itinerary of F^eis) 



Fi 



where N e 7 - 



= N < 31 



e (*+i)(j+i) 



ko fei 



(3) 



Consider periodic or- 



bit 00111 of a suspension of the standard horseshoe map 
equipped with the usual symbolic coding 2] (Figs. ^jp and 
Et). We find that as points gradually move in the section 
plane from their initial location to that of their image un- 
der the return map, triangle inversions occur when point 
4 successively crosses the four edges eis, ei 3 , e25 and 
e23. Thus the induced return map for edge itineraries is 
F 1 = Na^ 3 a^ 5 af 3 af 5 . For example, 



ei4e45 



N 



ei 4 e 45 — ► e 2 5e 5 i = Fi(ei 5 ), 



while edges not crossed by point 4 are trivially modified 
(e.g., ei4 625). This leads to the closed rule set 



ei4 — » e25, eis 



e25 esi, e 2 5 — > e 35 e 5 i, e 35 — ► e 4i (4) 



for edges in the invariant set of .Fi. Table [l] displays iter- 
ates F 1 m (ei 5 ) computed using |@}. Their length Fi m (w)| 
diverges exponentially as m — ► oo, indicating that trajec- 
tories in the neighborhood of the orbit are continuously 
stretched apart by the flow. The growth rate: 



HP) 



lim 

71— >00 



hi \F^(w)\ 



(5) 



is obtained as the logarithm of the leading eigenvalue 
of the transition matrix (M ee >), whose entries count oc- 
curences of edge e' or of its reverse in Fi(e) given by 0}. 
Here, /i(00111) ~ 0.5435. Table[Qalso shows that F* p (w) 
(p is the orbit period) converges to an infinite word 
satisfying Ff(w 00 ) — Woo, which is the analog of the in- 
finitely folded unstable manifold of the periodic orbit. 

The growth rate h(P) is expected to be the entropy 
1it{P) of orbit P, defined as the minimal topological en- 
tropy of a map containing P ||. Indeed, a piecewise 
linear map containing P with (M ee <) as Markov transi- 
tion matrix can be constructed and has entropy h(P), 
thus h T (P) < h(P). Conversely h(P) < h T {P), as h(P) 
is the minimal growth rate of the geometric length of 
curves passing through periodic points Pi and cannot be 
larger than the topological entropy of a map containing 
P, which is the supremum of stretching rates over curves 
in the plane . 

For a typical orbit, unlike in J3J, there are paths in 
the Fi-invariant set that trigger a "squeezing" rule 12bL 



(15) 

1 (251) 

2 (35152) 

3 (41525153) 

4 (5251535152514) 

5 (1535152514152515351525) 

6 (2514152515351525251535152514152515351) 

10 (1535152514152515351525251535152525153515 . . .) 

15 (1535152514152515351525251535152525153515 . . .) 

100 (1535152514152515351525251535152525153515 . . .) 

TABLE I: A few iterates F{™(eis) are given by their itinerary 
between periodic points [e.g., (35152) denotes the path 
e3ge5ieises2]. 




FIG. 5: (a) Periodic points of the horseshoe orbit 00111 and 
their schematic trajectory in section plane. Bold lines indicate 
edges involved in (b) Path P4P1P5P2P5P1P5P3 folds onto 
itself under action of induced return map Fi. The unimodal 
map obtained has 00111 as a periodic orbit. 



as for example eieee^ — » en for horseshoe orbit 0010111. 
Then Fi^ise^) ^ Fi(ei6)Fi(e67) and the transition ma- 
trix cannot be used for entropy computations, although 
estimates can still be obtained by direct iteration. In all 
examples we considered, enlarging the alphabet by recod- 
ing contracting paths as basis edges (e.g., ei67 = e\%e.&i) 
and applying other recodings required for consistency al- 
lowed us to rewrite Fi as an ordinary substitution like 10}. 
For example, the induced return map for horseshoe orbit 
0010111 can be rewritten as (e^fc = &ij£jk)'- 

ei4 e 2 5, eis — * 6257 ^76, &n — > ^257 e7i, e25 — > e 3 7 e^, 
637 — ► e 4 i, e67 — > e7i, ei67 — » e 2 5 e&i, e 2 57 — > e 3 7 e76i 

Besides ei67, basis path 6257 was introduced because its 
image overlaps ei67- A transition matrix can then be 
obtained, with entropy /i(0010111) ~ 0.4768. 

For all 746 periodic orbits of the horseshoe map up 
to period 12, we have compared growth rate JS} with 
topological entropy obtained by the train-track algo- 
rithm 0, E3, El- As illustrated in Table HU we found 
agreement to machine precision in each instance. This 
strongly suggests that in 3D, our approach is equivalent 
to the train track approach. Qualitative properties of 
chaos are also reproduced: the dynamics is determinis- 
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Orbit 



This work TTA 



01101 ?1 
001011 ?1 
00101 ?1 
001010 ?1 
001 ?1 

001110 ?1 
00111 51 

001111 °1 
001101 ?1 



0.4421 
0.3460 
0.4768 
0.4980 
0.5435 
0.4980 
0.4768 
0.3460 
0.4980 



0.4421 
0.3460 
0.4768 
0.4980 
0.5435 
0.4980 
0.4768 
0.3460 
0.4980 



Orbit 



This work TTA 



00010 ?1 
000101 ?1 
0001 ?1 
000111 ?1 

00011 ?1 

000010 ?1 
00001 ?1 

000011 ?1 
000001 ?1 



0.3822 
0.5686 
0.6329 
0.5686 
0.3822 
0.4589 
0.6662 
0.4589 
0.6804 



0.3822 
0.5686 
0.6329 
0.5686 
0.3822 
0.4589 
0.6662 
0.4589 
0.6804 



TABLE II: Topological entropies of positive-entropy horse- 
shoe orbits up to period 8 obtained with the approach de- 
scribed here and with the train-track algorithm (TTA). 



tic (by construction), invertible, and the stretching and 
squeezing processes are described in a symmetrical way. 

Remarkably, we note that while transformations © 
are invertible, the asymptotic dynamics is singular. Con- 
sider the itinerary wq = F 3 (ei5) = (41525153) in Ta- 
ble [I] which is the shortest subpath of visiting 
the four edges in As Fig. ]Ej[b) shows, the image 

Fi(w ) = (5251535152514)= (525153) + (35152514) con- 
sists of a subpath of wq concatenated with a reverse 
copy of Wo: this path is folded onto itself by a singular 
one- dimensional map. The same property holds for all 
subsequent iterates _F m (ei5), hence for the infinite word 
u>oo . This reflects that associated to an invertible return 
map (e.g., Henon map), there exists an underlying lower- 
dimensional noninvertible map (e.g., logistic map) de- 
scribing the dynamics along the unstable manifold, a key- 
stone of the Birman- Williams construction 0, . Note 
that the symbolic name 00111 can be recovered directly 
from Fig. [Sfb) using the usual coding for orbits of ID 
maps. This makes the new formalism promising for us- 



ing topological information to construct global symbolic 
codings as in 01 ■ How segments along wq are folded 
over each other and how neighboring cells are squeezed 
provide us with a combinatorial description of stretching 
and folding that could be used to determine the simplest 
template carrying the periodic orbit studied. 

To conclude, we have proposed that orientation preser- 
vation is a more general formulation of determinism than 
non-intersection of trajectories. In three dimensions we 
find that enforcing it on a triangulation of periodic points 
induces a nontrivial dynamics on paths along periodic 
periodic points. More precisely, a path map Fi is con- 
structed by: (i) following triangles advected by the flow 
as one rotates around the attractor, (ii) restoring orien- 
tation at each triangle inversion by exchanging opposing 
sides via transformations (J2J). When paths in the F\- 
invariant set do not experience contraction, entropy is 
obtained from an transition matrix indicating how ele- 
mentary edges in the invariant set are mapped among 
themselves. Otherwise, new basis paths must be intro- 
duced to account for contraction. A promising result is 
that despite its simplicity this formalism numerically pre- 
dicts the correct entropies for periodic orbits of the horse- 
shoe map. Preliminary calculations also suggest that it 
leads to a combinatorial description of the folding of the 
invariant unstable manifold over itself, yielding informa- 
tion about the symbolic dynamics of the orbit. It now 
remains to prove the validity of the approcach in 3D and 
to try to extend it to higher dimensions. 

This work grew out of innumerable discussions with 
R. Gilmore. I thank T. Hall, J. Los and F. Gautero for 
helpful explanations about train tracks, and M. Nizette, 
T. Tsankov, J.-C. Garreau, C. Szwaj and S. Bielawski 
for a careful reading of this manuscript. CERLA is sup- 
ported by the Ministere charge de la Recherche, Region 
Nord-Pas de Calais and FEDER. 
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